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We show how many-body ground state entanglement information may be extracted from sub- 
system energy measurements at zero temperature. A precise relation between entanglement and en- 
ergy fluctuations is demonstrated in the weak coupling limit. Examples are given with the two-state 
system and the harmonic oscillator, and energy probability distributions are calculated. Compar- 
isons made with recent qubit experiments show this type of measurement provides another method 
to quantify entanglement with the environment. 
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A standard assumption in thermodynamics is that the coupHng energy of the system to the thermodynamic bath 
in ' must be smaller than any other energy scale in the problem. In this paper, we explore the consequences of the 
. violation of this assumption when the combined system and bath are together in the overall ground state (or at zero 
temperature).^ From the thermodynamic point of view, this is a boring situation because nothing can happen: the 
system and bath cannot exchange energy. However, from a quantum mechanical point of view, the non-vanishing of 
Ch ' the coupling energy can play an important role for mesoscopic systems (where the thermodynamic limit cannot be 
I ^ ' applied). Thermodynamic relations must be applied only to the entire systemi^ In fact, even though the system 
is at zero temperature, if a measurement of a sub-system Hamiltonian is made, it can be found in an excited state 
^ with a probability that depends on the coupling to its environment. This non-intuitive result is a purely quantum 
I . phenomenon: it is a consequence of entanglement of the sub-system with the environment. In fact, we demonstrate 
' that knowledge of the probability to find the system in an excited state can be used to determine the degree of 
, entanglement of the sub-system and bath. Consequently, simple systems with well known isolated quantum mechanical 
^-H ' properties (such as the two-state system and harmonic oscillator) become "entanglement-meters" . 

There is growing interest in ground-state entanglement from the condensed matter physics community. Theoretical 
works on ground state entanglement have addressed entropy scaling in harmonic networks,^ spin-spin entanglement in 
quantum spin chainsS. and quantum phase transitions.^ * Entanglement properties of the ground state are also essential 
. in the field of adiabatic quantum computingiS Recently, there has also been interest in the relationship between energy 
frustration and entanglement >iS It is also interesting to link other ground state properties of a variety of mesoscopic 
systems to the zero-temperature entanglement energetics. These properties include the persistent current of small 
mesoscopic ringsiiiSii^ or of doubly connected Cooper pair boxes^^ii^*i^ single Cooper pair boxes measured by a 
, dc-SQUIDfiiiiS and the occupation of resonant statesii^ Furthermore, the role of entanglement with an unmonitored 
^ ■ environment in the decoherence of scattering quantum particles has been considered for many-body quantum chaotic 
^ ' bathg^° and recently at zero temperaturei^ 

It has long been recognized that the ground state properties of mesoscopic systems are very interesting. In particular, 
^ ' a small metallic loop penetrated by an Aharonov-Bohm flux exhibits a persistent current if the temperature is so low 
• i-H , that the phase coherence length becomes larger than the circumference. It is therefore of interest to investigate the 
' persistent current in rings coupled to a bath^^^ The ground state of a model of a ring with a quantum dot coupled 
H ' capacitively to a resistor was examined by Cedraschi et alm^ and it was found that the persistent current decreases 
■ - - t with increasing coupling strength and at the same time that the persistent current is not sharp but fluctuates with a 
variance that increases with increasing coupling strength. To explain these results these authors already alluded to 
energy fluctuations. Such an explanation implies a close connection between energy fluctuations and persistent current 
fluctuations. Indeed in the work presented here we substantiate this relationship. A simple and transparent model 
in which energy fluctuations can be investigated is that of an oscillator coupled to a bath of harmonic oscillators. 
Nagaev and one of the authors'^ calculated the variance of the energy of the oscillator as a function of the coupling 
strength to the bath. In the work presented here, we analyze not only the variance but the entire distribution function 
of energy of the oscillator in its ground state, and show how these fluctuations originate from entanglement. 

We consider a general Hamiltonian H — Hg + He + He, that couples (c) the system (s) we are interested in 
to a quantum environment {E) such as a network of harmonic oscillators ji^™!?*^ The lowest energy separable state 
is \S) = |0)s|0)£;, where |0){s_£;i, are the lowest uncoupled energy state of both systems. However, if the system 
Hamiltonian and the total Hamiltonian do not commute (which is the generic situation), then 15) is not an energy 
eigenstate of the total Hamiltonian. Thus, there must be a lower energy eigenstate (|0)) of the total Hamiltonian 
which is by definition an entangled state. Because time evolution is governed by the full Hamiltonian, the ground 
state expectation of any operator with no explicit time dependence will have no time evolution, insuring that any 
measurement outcome is static in time. This situation is in contrast to the usual starting point of assuming that the 
initial state is a separable state and studying how it becomes entangled. The reduced density operator of the system 
is given by tracing out the environmental degrees of freedom, p = Tr£;|0)(0|. Assuming the full state of the whole 
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FIG. 1: A mesoscopic ring with in-line quantum dot and Aharonov-Bohm fiux $, coupled to an environment comprised of an 
infinite chain of LC-oscillators. This system exhibits a persistent current in its ground state, which is related to the entanglement 
energetics. After Ref. ITU 



system is pure, the reduced density matrix contains all accessible system information, including entanglement of the 
system with its environment. Because repeated measurements of Hg will give different energies as the sub-system is 
not in an energy eigenstate, we are interested in a complete description of the statistical energy fluctuations. These 
fluctuations may be described in two equivalent ways. The first way is to find the diagonal density matrix elements in 
the basis where Hs is diagonal. These elements represent the probability to measure a particular excited state of Hs- 
A second way is to find all energy cumulants. A cumulant of arbitrary order may be calculated from the sub-system 
energy generating function, Z{x) ~ (exp(— xi/s)) (as always, (O) = TrpO) so that the n*'' energy cumulant is given 
by 

These cumulants give information about the measured energy distribution around the average. 

Before proceeding to calculate these energy fluctuations, we ask a general question about entanglement. Given the 
energy distribution function (the diagonal matrix elements of the density matrix only), can anything be said in general 
about the purity or entropy of the state? Surprisingly, because we are given the additional information that we are 
at zero temperature, the answer is yes. If we ever measure the sub-system's energy and find an excited energy, then 
we know the state is entangled. Although this statement alone links energy fluctuations with entanglement, a further 
quantitative statement may be made in the weak coupling limit. The reason for this is the following: the assumptions 
exponentially suppress higher states, so to first order in the coupling constant, we can consider a two-state system 
where the density matrix has the form 
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For vanishing coupling constant a = 0, the first term is just the density matrix for the separable state. The linear 
dependence of p on a holds to first order for the model systems considered below and is the entanglement contribution. 
If one measures the diagonal elements of p, one obtains pdown = 1 — pa and Pup — pa as the probability to be measured 
in the ground or excited state (because a is small, there is only a small probability to find the sub-system in the upper 
state). If we now diagonalize p, the eigenvalues are Ai^2 = {1 —pa, pa} -I- 0{a^). To first order in a, the eigenvalues 
are the diagonal matrix elements, so we may (to a good approximation) write the purity or entropy in terms of these 
probabilities even if the energy difference remains unknown. 

The Qubit. Let us now first evaluate the energy fluctuations of a qubit, a two-state system. The most general (trace 
1) spin density matrix is p — (JL + {po^Ox + {(^y)(^y + {'^z)'^z) /'^- A simple measure of the entanglement is given by 
the purity, Trp^ = {1/2){1+X'^ -f- -I- Z^), where Xi = {ai). It is well known that {X, Y, Z) form coordinates in the 
Block sphere. Purity lies at the surface where X^ + Y'^ + — 1, whereas corruption lies deep in the middle. 

We take the system Hamiltonian~ to be Hg — {e/2)az + {A/2)ax- Introducing the frequency = \/e^ + A'^/h 
and using the identity e"*'?" '^ = /cos^ — in ■ as'm^ with (3 = hx^, and the unit vector fi chosen to give e~''^^= 
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FIG. 2: Energy probabilities p^p and Pdoiun for the spin-boson problem. With increasing coupling to the environment it it is 
more likely to measure the qubit in the excited state. For the symmetric case (e = 0), we use the Bethe ansatz solution, while 
for the general case, we use a perturbative solution which is only valid for large e or large a. After Ref. Q. 



(n^ — e/Ml, Tlx — A/?ir2), it is straightforward to show 
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The energy probability distribution may be easily found by Fourier transforming Eq. , or by tracing in the diagonal 
basis of the system Hamiltonian. The answer may be expressed with only the average energy, {Hg) = f 
as a sum of delta functions at the system energies ±.Ml/2 with weights of the diagonal density matrix elements, 
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Clearly, if the spin is isolated from the environment, (Hg) = —hQ/2 (the ground state energy), the probability weight 
to be in an excited state vanishes. This distribution may also be found from knowledge of the isolated eigenenergies, 
the fact that (Hg) = EjPjj, and that Trp ~ 1. This later argument may be extended to n-state systems given the 
first n — 1 moments of the Hamiltonian and the n eigenenergies. 

Connection with Real Qubits. The probability weights depend on the energy parameters e and A, and the expecta- 
tion values of the Pauli matrices. For real qubits produced in the lab, these will depend on the environment. Often, 
we can link the basic phenomena we have been describing to physical measurements other than energy. Consider, for 
example, a mesoscopic ring threaded by an Aharonov-Bohm flux $ shown in Fig. ^ The ring has an in-line quantum 
dot coupled to it with tunneling contacts, where the tunneling matrix elements t^^tn depend on the flux $. Interac- 
tions between the ring and dot are described with the capacitances Cl and C^. The dot-ring structure is capacitively 
coupled to an external impedance Zf,xt modeled by an infinite chain of LC-oscillators. This external impedance plays 
the role of the quantum environment. The equilibrium state of the dot-ring system supports a persistent current as 
a function of flux. The persistent current is related to the effective two-level system operators only, and in turn may 
be related to the probability to find the excited energy state (for the symmetric case of e = 0), 
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where Io{^) is the uncoupled value of the persistent current. This physical implementation gives a direct translation 
between the measured persistent current and the entanglement between ring and dot. Different discussions of the 
effect of a bath on persistent current should be classified as whether the system Hamiltonian commutes with the total 
Hamiltonian (see Ref. .13.,) or does not (see Refs. 11,12). Another physical system that shows similar physics is a 
quantum dot connected in series with a tunnel junction, metallic reservoir and quantum impedance represented by a 
bosonic environment 

A common model for environmental effects is given by coupling the two-state system to a series of harmonic 
oscillators, the spin-boson modeliiiSiSiS^ In Fig. |21 we have plotted the upper and lower occupation probabilities 
for the spin-boson model as a function of the coupling constant a. For the symmetric case (e = 0), we have used the 
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Bethe ansatz solutioiiii*^^, while for e finite, we have used the perturbative solution in A/wc which is valid only for 
larger a or eiii Thus the plot is cut off at a small a. A computational _approach calculating the expectation values of 
the Pauli matrices over the whole parameter range was given in Ref. |2^- The quantum dynamics of this system was 
studied in Ref. '2?. One simple measure of the bath type is the slope of the occupation probability in the vicinity of 
a = 0. 

Experiments are always carried out at finite temperature, and it is important to demonstrate that there exists a 
cross-over temperature to the quantum behavior discussed here. For an isolated system in thermal equilibrium, where 
the coupling energy plays no role, the (low-temperature) thermal occupation probability is pth = exp[— (i?2 — Ei)/kT]. 
In the weak coupling limit for the symmetric spin boson problem, the probability to measure the excited state scales 
as Pup ~ — Q;log(A/LJc)-^"2, Setting these factors equal and solving for T* yields 

Since T* scales as the inverse logarithm of the coupling constant, it is experimentally possible to reach a regime 
where thermal excitation is negligible. If one carefully calculates many-body low temperature corrections to the zero 
temperature results, one obtain corrections quadratic in temperaturCfSS 

As an order of magnitude estimate, we compare with the Cooper pair bojii^ii^ which is among the most environ- 
mentally isolated solid state qubitsi^. Fromi^ which found a Q ~ 10'*, we estimate the quantum probability for the 
box to be measured in the excited state as Pup ^ — 10~*, which is of same order or larger than the thermal 
probability, pth ~ 10^^. Experimentally, Pup and pth may be confused by fitting data with an effective temperature, 
Pth oc exp(— /?off-ffs)^ However, one may distinguish true thermal behavior from the effect described here because 
Pup and Pth depend differently on tunable system parameters such as A. In fact, /3cff is an entanglement measure. 
The behavior discussed here is closely related to the breakdown of the concept of local temperature discussed in Ref. 

The Harmonic Oscillator. We now consider the entanglement energetics of a harmonic oscillator, Hs = p^ /(2m) + 
(l/2)ma;^(?^. Since there are an infinite number of states, the problem is harder. To simplify our task, we assume a lin- 
ear coupling with an harmonic oscillator bath. This implies that the density matrix is Gaussian so that environmental 
information is contained in the second moments {q^) and 
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Expectation values of higher powers of Hg are non-trivial because q and p do not commute. The purity of the density 
matrix Eq. ((JJ) is 

Trp2 ^ / dqdq'(q\p\q'){q'\p\q) ^ -=^= ■ (8) 



The uncertainty relation, •\/ {q^){p'^) > ?i/2, guarantees that Trp^ < 1, with the inequality becoming sharp if the 
oscillator is isolated from the environment. As the environment causes greater deviation from the Planck scale limit, 
the state loses purity. 

The generating function Z may be calculated conveniently by tracing in the position basis and inserting a complete 
set of position states between the operators, 

Z{X)= I dqdq'{qW){q'\e-^"'\q). (9) 



The first object in Eq. @ is the density matrix in position representation, given by Eq. Q. The second object may 
be interpreted as the uncoupled position-space propagator of the harmonic oscillator from position q to q' in time 
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{q'\e-^"^\q) = 



27r?i sinh fiujx 

This interpretation is quite general and may be used to extend this analysis to other systems. We find 

^, s { „ sinhey ... ^ l + coshevl ^ 

Z(x)= 2i? ^ + 2A(cosh£x-l) + ^ -\ , (11) 
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where e = hco, 2E = mu;'^{q'^) + {p'^)/m and A = 



E is the average energy of the oscillator, while A > 1 

0, 



is a measure of satisfaction of the uncertainty principle. Eq. has a pleasing limit for the free particle lu 

Z{x)free = {l+x{p^)/my^ , 



(12) 



which is just the generating function for Wick contractions, (p^") = (2n — 1)!! ((p^))". Thus, in Eq. 1111) . the inverse 
square root generates the right combinatorial factors under differentiation, and the nontrivial x dependence accounts 
for the commutation relations between q and p. The first few harmonic oscillator energy cumulants may now be 
straightforwardly found via Eq. 
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After inserting the mean square values for an ohmic bath (see the discussion above eqs. (|21I22I) '). Eq. (|13|l is identical 
to the main result of Nagaev and one of the authors 

Alternatively, we now consider the diagonal matrix elements p„„. An analytical expression for the density matrix 
in the energy basis may be found by using the wavefunctions of the harmonic oscillator. 
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where 7 = yj muj/fi and H„{x) is the n*'* Hermite polynomial^ In the energy basis, the density matrix is given 
by Pnm = J dqdq'ip^{q){q\p\q')'ipmi<l')- The position space integrals may be done using two different copies of the 
generating function for the Hermite polynomials/^ 
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The diagonal elements may be found by equating equal powers of the generating variables. We first define the 
dimensionless variables x — 27^(5^), y — 2(p^)/(7^?i^), and D = 1 + x + y + xy. x and y are related to the major 
and minor axes of an uncertainty ellipse. The isolated harmonic oscillator (in it's ground state) obeys two important 
properties: minimum uncertainty (in position and momentum) and equipartition of energy between average kinetic 
and potential energies. The influence of the environment causes deviations from these ideal behaviors which may be 
accounted for by introducing two new parameters, a — {y — x)/D, b = (xy — 1)/I? with — 1 < a < 1 and < 6 < 1. 
The deviation from equipartition of energy is measured by a, while the deviation from the ideal uncertainty relation 
is measured by b. We find 



= ^j-{b^-a^y'^Pn 



b/\/b^ - a? 
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where P„[z] are the Legendre polynomials. The first few energy probabilities are given below (without the \J AjD 
prefactor). 
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If we try and choose x and y so as to violate the uncertainty principle, unphysical results appear as some of the 
probabilities exceed 1, or become negative. The probabilities (|18|) also reveal environmental information. For example, 
Pii/poo = b a-iid is thus only sensitive to the area of the state, while P22/P00 = a^/2 + 6^ depends on both the 
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FIG. 3: The probability to measure a harmonic oscillator in the ground and first three excited states as a function of x and y 
(see text). The line traces out the behavior of the ohmic bath as a function of the coupling in the under-damped range. After 
Ref. Q. 



uncertainty and energy asymmetry. Additionally, if we expand the first density matrix eigenvalueS*22i with respect to 
small deviations of x and y, we recover pn in agreement with our general argument. To complete the circle, we may 
make an "energy transform" on these probabilities. 



oo 
n=0 



where En = (n + l/2)?ia; are the uncoupled energy eigenvalues of the harmonic oscillator. If we now identify the 
new generating variable t — \J\P' — a? cxp(— and deviation variable z — hj \/b^ — a? , we may make use of the 
summation formula for the Legendre polynomials;^ 

oo 

Y^rPn[z]^{l-2zt + f}-"^ , (20) 

ri=0 

to recover (after some algebra) the energy generating function Eq. 

Although X and y have been treated as independent variables, the kind of environment the system is coupled to 
replaces these variables with two functions of the coupling constant. For example, with the ohmic bath^»2^ (in the 
under-damped limit), the variables are 

. X 1 / 2 a \ 

= „ 1 - -arctan— == , (21) 

Vl — a V ^ Vl — a / 

y[a) = [l-2a^)x(a) + —\n—, (22) 

where a is the coupling to the environment in units of the oscillator frequency and lOc is a high frequency cutoff. 
This bath information is shown in Fig. Owith ujc = lOuj. The trajectory of the line over the surface shows how the 
probabilities evolve as the coupling a is increased from to 1 . Other kinds of environments would trace out different 
contours on the probability surface. 

In conclusion, we have shown that projective measurements of the system Hamiltonian at zero temperature reveal 
entanglement properties of the many-body quantum mechanical ground state. Consequently, repeated experiments 
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on simple quantum systems give information about the nature of the environment, the strength of the couphng and 
entanglement. The larger the energy fluctuations, the greater the entanglement. There are several possibilities for 
experimental implementations. We have mentioned measurement of persistent currenfeiiiS*i^ as well as projecting 
on the system's energy eigenstates. Another measurement possibility is a zero temperature activation-like process^ 
where the dominant mechanism is not tunneling, but the same quantum effects of the environment which we have 
discussed here. 
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